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Abstract 

A model for tlie motion of ring-shaped DNA in a gel is introduced and 

studied by numerical simulations and a mean-field approximation. The ring 

Tb ■ 

r^ , motion is mediated by finger-shaped loops (hernias) that move in an amoeba- 

o, 

, O,. like fashion around the gel obstructions. This constitutes an extension of 



previous reptation tube treatments. It is shown that tension is essential for 
describing the dynamics in the presence of hernias. It is included in the 



^\j , model as long range interactions over stretched DNA regions. The mobility 

o, 

t^^ \ of ring-shaped DNA is found to saturate much as in the well-studied case of 

G^ 

"^ \ linear DNA. Experiments in polymer gels, however, show that the mobility 

drops exponentially with the DNA ring size. This is commonly attributed 
^ , to dangling-ends in the gel that can impale the ring. The predictions of the 

o 

present model are expected to apply to artificial 2D obstacle arrays (W.D. 
k> \ Volkmuth, R.H. Austin, Nature 358,600 (1992)) which have no dangling-ends. 

H : 

, 5^ , In the zero-field case an exact solution of the model steady-state is obtained, 

and quantities such as the average ring size are calculated. An approximate 
treatment of the ring dynamics is given, and the diffusion coefficient is derived. 
The model is also discussed in the context of spontaneous symmetry breaking 



in one dimension. 
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1. Introduction 

Gel-electrophoresis is a widely used technique for separating DNA fragments according 
to size lip. The separation resolution is hmited by a saturation of the mobihty at large 
DNA size. Separation of large DNA fragments has been made possible by pulsed-field gel- 
electrophoresis [0,^ . In view of the phenomenal successes of these techniques, an analytic 
approach to the basic underlying motion of the molecule through the gel is desirable. 

Most theoretical treatments fl-[To| of the motion of the DNA through the gel are based 



on the reptation concept [^. The DNA is pictured as moving through an impenetrable tube 
defined by the surrounding gel obstructions, with the motion mediated by a snake-like rep- 
tation of the polymer ends. Reptation has proven very successful in describing equilibrium 



dynamics of polymers in gels and melts. Simulations [|12| and experiments [0, however. 



have indicated that for sufficiently long chains undergoing electrophoresis, an alternative 
mechanism of motion is important: the formation of finger- like hernias or leaks through the 
reptation tube. These hernias (sometimes also called hair-pins, loops or kinks) constitute a 
protrusion of the DNA chain through the walls of the reptation tube in a doubled-up loop. 
Hernias have been included in some recent simulations of linear DNA fragments undergoing 
gel-electrophoresis |T^-|IB|. An additional important effect, that is often neglected in treat- 
ments inspired by equilibrium reptation theory, is tension transmitted along the DNA chain 



|T^. Under a driving electric field, strong tension forces can dramatically affect the polymer 



motion 16,17 



In linear DNA chains, both hernia motion and ordinary reptation of the chain ends are 
possible. In order to separate out and emphasize the effect of hernias, we consider here DNA 
in the shape of a ring (open-circular DNA [0). The DNA ring can move around the gel 
obstacles only by hernias, sending out fingers in an amoeba-like fashion. There have been 
no theoretical studies on gel-electrophoresis of open-circular DNA, despite the fact that in 
practical applications, ring shaped DNA (plasmids) is often analyzed by gel-electrophoresis 
and shows behavior different from that of linear DNA fragments 



The behavior of ring polymers in the absence of an electric field is also of interest |23|-|25 



This problem is related to the behavior of a melt of ring polymers, and also to electrophoresis 
in the weak-field limit through an Einstein relation. The diffusion of ring-shaped polymers in 
a lattice of obstructions has been treated by numerical simulations and theoretical arguments 



25| . An exact treatment of the statics and especially the dynamics of ring-shaped polymers 
in zero field is, however, not available. 

In this work, a model for the motion of ring-shaped DNA in a gel is introduced and 
studied numerically and analytically. The ring motion is mediated by finger-shaped loops 
(hernias) that move in an amoeba-like fashion around the gel obstructions. This model, 
described in Sec. 2, constitutes an extension of previous reptation tube treatments. It is 
instructive to first study the model neglecting the effects of tension transmitted along the 
DNA polymer. Monte-Carlo simulations of the model, summarized in Sec 3a, show that 
the chain mobility in this case decreases exponentially with DNA size. This is due to the 
formation of hooks which reduce the mobility. This behavior is modified when tension is 
taken into account. In Sec. 3b, tension is added to the model as long-range interactions 
over stretched regions of the chain. Tension increases the unhooking rates and stabilizes a 
ring conformation aligned with the field. This causes the mobility of long ring-shaped DNA 
to saturate to a finite value, much as in the well-studied case of linear DNA. Experiments in 
polymer gels [|T9| -^, however, show that the mobility drops to zero with the DNA ring size, 
with large rings hardly penetrating into the gel. This is commonly attributed to the rings 
becoming impaled on dangling-ends in the gel. The predictions of the present model are 



expected to apply to artificial 2D obstacle arrays p8[ which have no dangling-ends. In Sec 



3c, the polymer motion is qualitatively described by a mean-field treatment. In the zero- 
field case, discussed in Sec. 4, an exact solution of the model steady-state is obtained, and 
quantities such as the average ring size are calculated. In Sec 4b, an approximate treatment 
of the zero-field ring dynamics is given, and the diffusion coefficient is derived. This gives 
an analytic foundation to previous scaling arguments |^, and suggests a framework for 
analysis of dynamical features of driven polymers. In Sec. 5, the model is also discussed in 
the context of spontaneous symmetry breaking in one dimension. 



2. Model for DNA in a Gel Including Hernias 

We present a model for a charged polymer ring reptating in an electric field in an array 
of obstacles (eg. a gel). The model is based on the Rubinstein-Duke (RD) approach 10-3 , 
and is extended here to include hernias, which are hair-pin shaped excursions out of the 
usual reptation tube. Hernias are crucial for polymers in the shape of a ring, in which the 
motion through the surrounding obstacles may be accomplished only by hernia fingering. 

We begin by describing the RD model for reptating linear polymers. We then extend 
the model to include hernias. In the RD model, the gel is idealized as a lattice of point 
obstacles, with pore diameter b, as shown in Fig 1. In agarose gels, b ~ lOOnm, while in 
recently introduced artificial obstacle arrays b ~ Ifim [^. The DNA is represented as a 
chain of L segments, each of one persistence length (~ 50n?Ti). The segments may either be 
stretched, when the polymer threads through adjacent pores, or coiled, when the segment 
is contained in one pore. Each configuration of the polymer is represented by the positions 
of the successive cells that the polymer threads. A simplified description of the chain is 
coded by the projection along the field direction of the displacement between segment ends. 
This displacement can be either +b, when the segment threads between two pores in the 
field direction, —b when it threads two pores against the field direction, or when the 
segment is coiled in one pore. Thus the configuration is reduced to a ID lattice of L sites. 
Each site i corresponds to a DNA segment, and has a state (pi, which can be either -|-, — 
or 0, as demonstrated in Fig 1. Note that in this description some information regarding 
the microscopic configuration is lost. However, it provides a convenient way to model the 
dynamics [0-3 • 

The chain reptates by the motion of the coiled, lax segments through the chain. In 
aquaeos solution, the DNA is assumed to be uniformly charged. The forces acting on each 
segment are an electric force F^. = QE where Q is the charge per segment and E is the 
field strength, and a thermal Brownian noise Fth of the order of kT/b. These forces are 
represented in the model by the following rules. At each time step, a pair of sites is chosen 
at random, and a move is made with the following rates: 



+ — > 0+, at rate q, OH > +0 at rate p (1) 

—0^0— at rate p, > —0 at rate q (2) 

while H — or — h pairs are stuck, since they represent two stretched segments hooked around 
a gel obstacle (see Fig. 1). Moves in the field direction are favorably biased over the reverse 
moves, through the rates p and q. These rates are determined, in the case of weak fields, by 
local detailed-balance conditions f^, such as 

p = uJoexp{e/2) (3) 

q = uJoexp{-e/2) (4) 

where ljq is a microscopic rate and e = QEb/kT. Note that this is a nonequilibrium dynamics 
and it does not obey full detailed-balance. The ratio between these rates is thus a Boltzmann 
factor of the ratio of electrical to thermal energy, 

p/q = exp{e) e = QEb/kT (5) 

These rules, along with rates for injection of + and — at the head and tail of the linear 
chain, define the Rubinstein-Duke model [§|,0. 

We now extend this model to account for hernias. A hernia amounts to a projection 
of the chain from a pore with at least 2 coiled segments (2 adjacent sites) into a new 
pore, threading one segment out of the pore and another segment back to the original pore 
(for example move B in Fig. 1). This corresponds to a pair creation move 00 —>■ H — . 

The reverse annihilation move H > 00 corresponds to the hernia retracting and forming 

two coiled segments in a single pore. After a pair is created, the + and — can diffuse 
according to the RD rules. An important feature of the model is that pairs of + and — 
which are created together are tracked as a connected pair throughout the dynamics. Each 
-|- in the configuration has a unique — to which it is connected. Keeping track of such 
connections between pairs is necessary in order to track the hernia finger hierarchy. To 
see this, consider a pore with many coiled segments. A number of hernias may be formed, 

6 



projecting into different neighboring pores. An important point is that +'s and — 's from 
different hernias can not annihilate (assuming that several hernias originating from the a 
pore always project to different neighbors, a reasonable assumption for pore lattices of high 
coordination number). Thus, the pairings of +'s and — 's must be tracked: each + can 
annihilate only with the unique — to which it is paired. Each configuration is defined by the 
+ and — and sites, along with their pairing to hernias (Fig. 2). Only pairings in which 
the hernia pairs are nested are allowed, as shown in Fig. 3 (pairing lines may not cross 
each other). Starting with an allowed configuration, the hernia creation rules assure that 
the configuration at each subsequent time is also allowed. The phase space is larger than in 
the original RD model which has the three states +,— and for each site, but no pairings. 
The hernia creation and annihilation moves, which supplement the reptation moves of 
Eqs. m and 0, are: 

00 — > +u— , at rate c, 00 -^ — u + at rate d (6) 

+u > 00, at rate a, — u H > 00 at rate a! (7) 

The symbol u denotes pairs of + and — that have been created together (connected pair). 
The hernias tend to grow, as the field bias pushes the +'s to the left and — 's to the right. 
The hernias may develop sub-hernias, and a hierarchy of hernias may form. An example is 
shown in Fig. 2, for a ring shaped polymer. Thus, the polymer can assume a highly ramified 
shape, with a hierarchy of hernias of different sizes. 

The model as described above neglects an important physical effects: the tension trans- 
mitted along the chain. As shown below, this proves to be very important in the elec- 
trophoresis of ring shaped DNA. Tension acts as a long range effective interaction, and is 
included in the model as described in Sec 3b. 
3 Gel Electrophoresis of Open-Circular DNA 



Using the model, we studied gel-electrophoresis of ring-shaped (open circular [|T^) DNA. 
Periodic boundary conditions are thus imposed in the model. The ring is not concatenated 
with any gel obstacle (it is prepared outside the gel and moves into the gel under the field 



influence). We first study the model in the absence of tension by Monte-Carlo simulations. 
The treatment of tension, and its effects on the dynamics are presented in Sec 3b. A simple 
mean-field treatment is given in Sec 3c. 

3a. Monte-Carlo results: It is instructive to first study the model as described in Sec 2, 
neglecting the effects of tension. In order investigate the model, we performed Monte-Carlo 
simulations, typically using a = c = 1, a' = p,c' = 0, p/q = 1.01 — 2, and ring lengths up to 
L = 100. The mobility as a function of time is shown in Fig. 4. It is seen, that the mobility 
displays a spiked behavior: the system is effectively in one of two states, one with a positive 
mobility, and one with a zero mean mobility. The average lifetime of the zero-mobility states 
grows as L increases. 

The nature of the dynamics is clarified by snapshots of the ring configurations in the 
two states, shown in Fig. 5. It is seen that the ring cycles between quasi-linear and hooked 
states. The high-mobility phases corresponds to the quasi-linear conformation in which the 
ring is aligned with the field [^ . The conformation of effective charges in the model that 
represents this conformation is shown in Fig. 6a. This conformation is short lived, because 
it develops an instability: a hernia that buds on the side of the quasi-linear ring grows into 
a hooked configuration with two stretched arms, pinned over an obstacle (see Fig. 6b). In 
this phase, the ring is stuck, and there is zero mean mobility. The hooked state persists 
for a long time. The ring unhooks by one of the arms retracting by fluctuations, until a 
new, quasi-linear, high mobility shape is attained. This explains the burst-like structure of 
the mobility. In the absence of tension, the mean mobility /i decreases exponentially with 
L (Fig. 7). This is because the unhooking rate, by which an unstretched segment moves 
from one arm to the other, is exponentially small, since it takes of order L steps against 
the field for the segment to escape the arm's effective potential trap. Similar dynamics can 
occur also in linear chains [|1^,|T^. In the unhooking phase, tension plays a crucial role, as 



described in the following section. 

3b. Effect of Chain Tension: It is important to consider the effects of tension 
transmitted along the chain |ll2| , p!6| , p!7| , |30|l . The main effect of tension is to dramatically 



8 



increase the unhooking rates of stretched hooks. It acts as a long-range interaction between 
coiled segments |T6[ . Tension in the context of gel-electrophoresis of linear DNA was treated 
in a previous study |[T^, where coiled segments were allowed to make long-ranged hops along 
the chain. The present treatment simplifies this by using only local hops. In addition, the 
present model extends Ref []TB[ by taking into account the effect of tension on the hernia 



creation and annihilation rates. 

In order to model the effect of tension, we note that under the influence of an electric 
field the charged chain behaves like a chain moving in a gravitational field coupled to its 
weight. The tension transmitted along the chain is relaxed at coiled (0) segments and at 
hernia tips (+u— and — u+ paired at neighboring sites). Note that unpaired neighboring H — 
or — h sites represent segments of DNA chain which are draped over a gel obstruction and 
therefore are capable of transmitting tension. We define an effective field for each pair of 
sites, which corresponds to the tension generated by regions of stretched chain on both sides 
of the sites. This results in a movement rate Uj for the pair of sites j and j + 1, which depends 
on long-range interaction between different sites, as described below. At each step of the 
simulation, a pair of sites j and j + 1 is chosen at random. Pairs at which there is an allowed 
move are of three types: (a) coiled site adjacent to a stretched site (f)j,(f)j^i = 0—, —0,0+ 
or +0, (b) two coiled sites 4>j, 0j+i = 00 (creation move) or (c) a hernia tip 0j, 0j+i = +u — 
or — u+ (annihilation move). We will refer to such pairs as relaxed pairs. If the pair is not 
one of these three types, it remains unchanged. If the pair has an allowed move, the move 
is performed with the rate Wj, and the states of sites j,j + 1 are accordingly adjusted. A 
new pair is chosen and the process is repeated. 

To derive the movement rate for pair j (site j and j + 1), Uj, we consider the tension 
transmitted along the DNA due to stretched regions of chain on either side of the pair. 
Since the tension accumulates along these regions, the local tension field is proportional to 
the net displacement in the field direction of these stretched regions. The stretched regions 
terminate at either a coiled segment (0 site) or a hernia tip since theses are points when the 
chain tension is relaxed. The effective tension force acting on a pair consisting of a coiled 
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segment adjacent to a stretched one (O,0j+i), with 0^+1 = ±1, is 

F3 = \^ E <t^m (8) 

•m=j+l 

where ki is the closest succeding site to site j + 1 which is a member of a relaxed pair, and 
the dimensionless external field is e = QEb/kT. Similarly, for a (0j,O) pair with (pj = ±1, 
the effective force is given by 

Fj = -^e J2 <Pm (9) 

m=k2 

where /c2 is the closest preceding site to site j which is a member of a relaxed pair. The 
force acting on a hernia tip (+u~ or — u+ pairs) is 

F, = l<fl K-Y. <^rn) (10) 

m=j+l m=k2 

with similar definitions of ki and ^2. (We note that the model can also be applied to linear 
chains, where additional points at which tension is relaxed are the chain ends). 

As an example, consider the configuration of Fig 2. This configuration contains one 
coiled segment at site i = 10, 0io = 0. Consider sites 9 and 10, at which the configuration is 
—0. To evaluate the effective field for this pair, we sum in the two sites to the right of the 
pair (these sites cancel each other), where we reach a hernia tip. The total effective force 
is Fg = 0. At sites 10 and 11, where the configuration is 0+, tension accumulates along a 
three sites stretched + region to the right of the pair, which terminates at a hernia tip, and 
Fio = |e. 

The movement rate ujj can be related to the local tension force Fj from a consideration of 
the thermal and friction forces on the string. The motion of the DNA segments through the 
solvent is such that viscous drag forces are much larger than any inertial term [^,|1^]. The 
behavior of the chain under thermal noise can be treated using a Fokker-Planck approach 
||16|| , using a Smoluchowski equation for the string motion along the tube contour, which 
includes a Brownian term and a friction coefficient proportional to the string's length. In 
the present work we propose a simpler physical model, which is valid at both the strong and 
weak field limits: 
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,exp(F,) F,<0 
u;,= (11) 

r,(l + F,) F,>0 

where Tj is equal to a microscopic rate cuq for pairs where a coiled segment can move 

(0j, 0J+1 = +0, 0+, 0—, —0), Tj = Co for pairs where a hernia can be created ( (pj, 0j+i = 00) 

and Vj = tto where a hernia can be annihilated {(j)j,(f)j^i = +u~; ~u+ )• At all other 

pairs of sites, Tj = 0, since no other moves are allowed. The constants ujq,Cq and oq are 

the microscopic rates of the various processes. Eq. |TT| goes to a Boltzmann factor for low 

effective fields where it represents local detailed-balance. At high effective field strength, 

the movement rate becomes linear in the effective field strength. This is expected since at 

large fields thermal fluctuations become unimportant and the local chain velocity becomes 



proportional to the local force [|l^]. Eq. |TI] allows hernia annihilation rates to be affected 
by chain tension, with annihilations at hernia tips flanked by long regions stretched in the 
field direction given a high rate. 

The model without tension, described in Sec. 2, can be recovered from this model, by 
taking ki = j , k2 = j + I in Eq. |8|-|T0|. This corresponds to the screening limit, when 



the density of coiled segments is so high that there appear no extended regions of stretched 
segments in which tension can develop, and the field at each bond is due to the external 
field alone. In this case, the movement rates are related to those of the model without 
tension described in Sec 2 at low fields (e << 1) via: p = uj^e"!"^ ^ q = cuoc""/^, c = cq, 
a = aoe~'^,c' = co,a' = a^e'^. It is seen that in this limit, the rates satisfy Eq. |[ 

Tension causes hooks to have a much smaller effect on the mobility. Monte-Carlo calcu- 
lation including tension at two different field strengths are shown in Fig. 7. We find that the 
ring mobility saturates at large DNA sizes, much as in the well studied case of linear DNA 
fragments. It is interesting to note that the ring arranges itself into a quasi-linear shape in 
the size regime studied. Here, coiled segments are very frequent - roughly 1/3 of the seg- 
ments are coiled. Thus, tension is screened by the coiled segments and has a very small effect 
during most of the dynamics, since it is important only in long, continually stretched pieces 
of the chain. Only when a side-hernia forms, and a hook begins to be created, does tension 
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come into play, and essentially stabilizes the quasi-linear shape aligned with the field. We 
note that for very large rings, branching effects similar to those discovered in large linear 



fragments in Ref |T6|, are likely to occur, though the ring mobility should remain constant. 

The screening of tension by coiled segments is very important in explaining the sucsses 

of reptation tube models for linear DNA which seems to describe experiments on linear 



DNA fragments quite well |]T0|,^, though the models neglect both hernias and tension. The 
present results suggest that when including hernias in a model of polymer dynamics, it is 
essential to also take tension into account, as these two effects have a canceling behavior, 
respectively increasing and decreasing the hooking rates. 

The predictions of the present model that the mobility saturates with the ring size seem 
to disagree with experiments. Studies of open-circular DNA (plasmids) run through agarose 
gels show that above a certain DNA size, the plasmids are "stuck at the wells" and do not 



enter the gel ||T9|-p2|. The explanation offered by Micel et. al. |T9| is that the rings become 
hooked on dangling ends in the gel (unconnected ends of the gel fibers or other impurities 
that penetrate the pores), which impale the ring ("hoop in stick" effect). The ring can 
unhook by a fluctuation which can overcome the fleld pulling the ring. The probability of 
such a fluctuation is exponentially small in the ratio between the electric force pulling the 
ring and the thermal forces, and the mobility is 

/i ~ exp{-QENb/kT) (12) 

The saturation of the mobility predicted in the present model could be checked experimen- 



tally on recently introduced artiflcial 2D arrays of obstacles [^ with no dangling-ends that 
can impale the ring, as suggested in Sec. 6. 

3c. Mean- field treatment: 

In the presence of tension, the DNA is found mostly in a quasi-linear conformation 
aligned with the fleld, with many coiled segments. The coiled segments essentially screen 
tension. This allows for a simple and local mean-fleld treatment of the DNA motion. 

Consider a quasi-linear chain, (Fig. 5, rightmost and leftmost conflgurations) . In this 
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configuration, hernias are annihilated at the upheld end of the ring. The coiled segments 
(O's) which are generated move down to the leading end, where a new hernia is formed. 
The density of coiled segments, p, is given by a balance of hernia creation and annihilation. 
The rates for these processes are obtained in the mean field approximation by neglecting 
correlations: Hernias are created (at the leading end) when two coiled segments are adjacent, 
at a rate cp^, and annihilated (at the upfield end) when two stretched segments are adjacent 
at rate a(l — p)^. In this approximation for the annihilation rate, the assumption that the 
ring is quasi-linear is used, since a stretched segment can annihilate only with its unique 
pair. In a random configuration the pair would have a small chance of being adjacent. Here 
we assume that in the quasi-linear configuration, a pair of stretched segments at the leading 
end may always be annihilated. 

The balance between annihilation and creation yields 

p = 1/(1 + ^a) (13) 

The mobility /i, equal to the mean center of mass displacement per unit time, is given by 
an average over all the configurations allowing movement, weighted by the respective rates. 
Since in steady state the hernia annihilation and creation moves balance each other, we have, 
in the simplest mean field approximation that the mobility is proportional to the probability 
of finding a stretched segment adjacent to a coiled one: 

p={p-q)p{l-p) p = l/(l + V^) (14) 

The qualitative features of the mobility are reasonably described by the simple mean- 
field theory as shown in Fig. 8, where the density of coiled segments and the mobility as 
a function of the ratio of hernia creation and annihilation rates c/a are shown. At high 
ratio of creation to annihilation rates c/a, the chain is dense with stretched segments and 
the mobility is low. At low c/a, there are few stretched segments that can move, and the 
mobility is again low. Around c/a = 1, where the density of coiled segments is around 1/2, 
the mobility is at a peak. The simulations show similar qualitative behavior. The mean-field 
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mobility overestimates the full model mobility by a factor of about 2. This is probably due 
to processes that impede the motion, such as hooking and pair creation in the bulk of the 
chain and not only at the head and tail, that are neglected in the mean-field treatment. 
4. The Zero Field Case 

We now turn to the case of zero electric field. This case is important as a question in 
classical polymer physics p3| , ^ : what is the effect of hernias on the statics and dynamics 
of a chain in a gel or melt at equilibrium? In addition, the zero-field diffusion can be related 
to the low-field electrophoretic mobility via the Einstein-Nernst relations. 

The zero-field case offers a significant simplification in the model: tension can be ignored 
in this case, and the model described in Sec 2 is used, with q = p, c = c' and a = a'. The 
probability of a given configuration C, P{C), is governed by the Master equation 

dP{C)/dt = J2{^iC' -^ C)P{C') - W[C -^ C')P{C)} (15) 

c 

where W(^A —* B) is the rate of transition from configuration A to B. A solution to the 

master equation is found which satisfies detailed-balance. Each move which preserves the 

number of stretched segments, such as +0 -^ 0+, is exactly balanced by its reverse move. 

Moves where hernias are created or annihilated are balanced by the reverse move, with an 

extra factor related to the creation and annihilation rates. A configuration A which has a 

00 at a certain bond, can, in a single move, go to or be reached from only two configurations 

B, B' which are exactly the same as A except that they have either a +u— or a — u+ at the 

bond. The solution for the probability of configuration C is 

P{C)=N{L)-\c/a)''^^^ (16) 

where c/a is the ratio between hernia creation and annihilation rates, and h{C) is the total 
number of hernias in the configuration C. In the steady state, all configurations have equal 
probability, up to a factor depending only on the total number of hernias in the configuration. 
This solution is remarkable in that although there are strong interactions between different 
hernias, the probability of each configuration depends only on the number of hernias and 
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not their relative positions and sizes. The normahzation factor N{L) = J2ci'^/'^)'^^'"^ i^ 
connected to the number of allowed configurations (only configurations with nested hernias 
are allowed, as shown in Fig. 3). N{L) satisfies the recursion relation: 

L-2 

N{L) = N{L - 1) + {2c/ a) ^ N{l)N{L -2-1) (17) 

The terms on the right-hand side can be understood as follows: given a ring of size L, 
choose a site. The first term corresponds to the case where the site contains a 0, and 
thus the configurations of hernias can be mapped to a ring of size L — 1 by deleting the 
site. The second term is the case where the site is a member of a hernia with its pair 
at a distance of / + 1 sites (hernia of size /). The factor 2 in Eq. ^^ is due to the two 
possible assignments of + and — charges to a hernia pair, which have equal probability in 
the absence of an electric field. The values iV(0) = iV(l) = 1 are supplemented to this 
recursion relation. Solving for the asymptotic form for N{L) at L ^ 1, we consider the 
function u{L) = ri~^N{L). For sufficiently large t] it has a finite integral. Using a Laplace 
transform of u{L), g{s) = J2'L=o'^i^)^~^^^ ^^ the recursion relation Eq. [1^, we find p^ 



1 - e-'r]-^ + J{e-'r]-^ - 1)2 - 4:{2c/a)e-^-'T]~^ 
'^'^ = 2(2c/a),-2e-2- ^'^^ 



The smallest value of rj for which g{0) = Y, u{L) exists is rj = 2J2c/a + 1. At this value of 
r], for small s, g{s) '^ gQ + gis^^"^- This corresponds to the following asymptotic form of the 
partition sum at L ^ 1 



N{L) = NoL-^/\2^2c/a + if (19) 

with 



No = {4^c/a)-\l + 2^2c/a)[^2c/a + Ac/a]^/^ (20) 

This allows derivation of steady state densities, such as p, the density of coiled segments 
(O's). This density is given from the construction of the recursion relation Eq |T^ simply 
by those configurations at which a given site is 0, compared to the total weight of the 
configurations: 
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p = N{L- l)/N{L) = (2^2^ + 1)-^ (21) 

The hernia size distribution, n{l), defined as the probabihty that a selected site belongs to 
a hernia of size /, is n{l) = N{l)N{L - 2 - l)/N{L) (see Eq O). Thus, for 1< / < L/2, 

n{l) ^ NqI-'^'^ (22) 

This suggests that the ring polymer adopts a ramified fingered shape, with a power-law 
spectrum of finger sizes. 

4a. Average Ring Size: The exact solution also allows calculation of the mean size of 
the ring at zero field. The radius of gyration serves as a convenient measure of the ring size: 
it is defined as the mean squared distance between two sites (i.e. the sum of number of + 
sites minus the number of — sites between i and j squared, averaged over all configurations 
and all i and j). To evaluate the radius, we define the function 

G{i) = Y.RU^P{C) (23) 

c 

where 

RUi) = {J:^kf (24) 

fc=l 

where (pk is the state of site k in configuration C . That is, G{i) is the average over all 

configurations of the squared distance between the two "test sites" 1 and i. The radius of 

gyration R is simply given by the root of the mean of G{i): 

R' = L-'J:G{z) (25) 

1=1 

We construct a recursion relation for G{i), noting that the only contribution that does not 
average to zero is from unclosed hernias in the interval between the sites. Since there is no 
field, each site on a hernia can be a + or a — with equal probability. The mean squared 
distance, averaged over all assignments of + and — to the hernias, is just the number of 
unclosed hernias. To build a recursion relation for G{i), we go from an interval of size i — 1 
to size i. Thus, to the interval of sites l...i — 1 we add one site, i, which may be either a 0, 
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belong to a hernia that closes outside the interval, or close one of the open hernias in the 
interval. Thus, we divide all of the possible configurations into three groups: Cq in which 
site i is in a state, Ci in which site i belongs to a hernia that closes outside the interval, and 
configurations C2 in which site i pairs with a site inside the interval, thus closing a hernia. 
In configurations Cq, R%{i) = Rci"^ ~ !)■ ^^ configurations Ci, R%{i) = R^ii — 1) + 1. In 
configurations C2, R^{i) = R^ii — 1) — 1. Thus, 

G{i) = Y: RU^ - l)^(C) + E PiC) - E PiC) (26) 

C Ci C2 

This leads to the recursion relation 

L-2 

G(i) = G{i - 1) + (2c/a)[ J2 N{l)N{L -2-1) (27) 

l=i-l 

-Y.m)N{L-2-l)]/N{L) 
1=0 

The boundary condition is G{1) = 1 — N{L — 1)/N{L), since the mean squared displacement 
due to a single site is just the density of uncoiled segments 1 — p. G{i) rises to a peak when 
i = L/2 and then drops off to at 2 = L, because the ring is closed and going around the 
ring the mean displacement must go to zero. Going to the continuum limit, and using the 
symmetry of G{i) around i = L/2, we find 

dG{i)/di = {2c/a)N{L)'^ f " N{l)N{L - l)dl (28) 

For large z, using the asymptotic form of N{L) found above (Eq. |T9|), the integrals involved 
in calculating G{i) can be performed analytically. This yields 

G{i) = 8{2c/a)NoViVL^/VL (29) 

showing that the mean maximal excursion of the ring is G{L/2) oc -\/L- Integrating over 
G{i), one obtains 

R = ^■K{2c/a)NQL^'^ (30) 

This scaling is valid for long chains. The radius for finite chains can be readily found from the 



recursion relations |1^ and ^. Thus, the ring adopts a much more compact configuration than 
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the linear, reptating chain, in which R oc L^/^. This form is vahd for chains small enough 



that excluded volume effects can be neglected |25,31|. An intuitive argument predicted the 



R ~ L^/^ scaling [^,^, by mapping the ring to a randomly branched graph. The exact 
solution of present model allowed us to derive this scaling exactly, from a steady state of a 
dynamic model. 

4b. Dynamics of a Ring: The solution for the steady state allows a rather accurate 
approximation for the polymer dynamics. We consider the motion of a marked hernia on 
the ring. We define P{l,t) as the probability of finding the hernia at size / at time t. The 
hernia can grow if there is a adjacent to the hernia from the outside, and it can shrink if 
there is a zero adjacent to the hernia from inside. We make the approximation that these 
probabilities are given by the corresponding steady-state probabilities. This approximation 
is a good one for large hernias. This is because the motion of the hernias is diffusive, and 
hence many configurations of the ring segments inside and outside the marked hernia are 
sampled on the time-scale of the effective hernia motion. The evolution equation for the size 
of the marked hernia P{1, t) is 

dP{l, t)/dt = Pi{l + 1)P{1 + 1, t) + Poil - l)P{l - 1, t) (31) 

-[A(/)+p,(/)]P(/,t) 

where pj(/), Po{l) are the probability of a zero adjacent to the size-1 hernia from the inside 
and outside respectively (a microscopic rate constant has been factored out of the equation 
so that each move take one time unit). At the boundaries, the equations have terms cor- 
responding to the annihilation of the hernia, which is possible when it is of size / = or 
l = L-2: 

dPiO,t)/dt = A(l)P(l,t) - [a + po(0)]P(0,t) (32) 

dP{L - 2, t)/dt = po{L - 3)P{L - 3, t) (33) 

_[a + p.(L-2)]P(L-2,t) 

The main idea of the present approximation is to use the exact steady state solution to 



estimate the probabilities of motion for the hernia. This yields 

p,{l) = N{1 - 1)/N{1) ~ p(l + h-') (34) 

Po{l) = N{L-l- 1)/N{L - /) ~ p[l + ^{L - I)-'] (35) 



where the asymptotic forms are valid at / ^ 1 and L — l^l, and p = (2w2c/a+ 1)~^. Note 
that smaller hernias have on average more zeros inside them than larger hernias. This is an 
entropic effect, due to the larger number of hernia-pairings in a large hernia. This creates 
a bias for small hernias to shrink, an effect that has important consequences for the ring 
dynamics, as shown below. 

In order to analyze the scaling of Eq. |3^, it is convenient to turn to a continuum form: 

dP{l, t)/dt = pV^P{l, t) + pV[f/(/)P(/, t)] (36) 

supplemented with appropriate boundary conditions (the exact form of the boundary con- 
ditions do not affect the scaling results obtained below). The potential 

U{l) = ^-[l-'-{L-ir^] (37) 

corresponds to the bias of small hernias to shrink. Consider the case of a marked hernia 
created at time t = 0. Thus, the initial conditions are P{l,t = 0) = 6{l = 0). We define 

/■oo 

F{1) = / P{l,t)dt (38) 

Jo 

/■oo 

r(/) = / tp{i,t)dt (39) 

Thus the mean lifetime of the marked hernia is 

/■OO /"OO 

r = / T{l)dl/ / F{l)dl (40) 

Jo Jq 

Ordinary differential equations for T{1) and F{1) may be easily formed by appropriately 
integrating over Eq ^: 

- p(/, t = 0) = V^F{1) + V[f/(/)F(/)] (41) 

-F{1) = V^T{1) + V[f/(/)T(/)] (42) 
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These equations are exactly soluble, yielding rather complicated expressions. However, to 
understand the scaling form of the solutions, it is useful to consider a simpler potential of 
the form U = u^/l, which is equivalent to the full potential Eq. ^at / << L, with uq = 3/2. 
Plugging in power law forms for F and T in Eq. |¥T| , ^ , we find F{1) ~ l~^° and T ~ l'^~'^°. 
For Mo = 3/2, /o°° F{l)dl ~ F(0) and J^T{l)dl ~ L^/\ yielding for long chains 

r ~ roL^/2 (43) 

The same asymptotic result is found using the full potential of Eq. ^ The lifetime of a 
given marked hernia is much shorter than that expected from a simple diffusion argument, 
T r^ L"^ (which corresponds to using no potential, U = 0, in the above calculation). This is 
due to the entropic bias of small hernias to shrink, which compels hernias to spend less time 
in "random walk" motion as compared to the pure diffusion case. Note that the power law 
exponent 3/2 in Eq. ^ is related directly to the prefactor of the effective potential that is 
derived from the detailed steady state solution. 

From this result, one can readily derive the scaling of the ring center of mass diffusion 
coefficient, D, using a classical scaling argument ||TT|. Essentially, diffusion proceeds by the 
transport of hernias along the ring. Each hernia takes a time r ~ L^ to travel a distance of 
order R{L) ~ L'^, the linear size of the polymer. If there were only one hernia on the chain 
at any time, the center of mass diffusion constant Dq would obey 

Dor = {R{L)/Lf (44) 

where the right hand side represents the mean square displacement of the center of mass 
arising from the transport of one hernia across a distance R. In reality, the number of hernias 
present at any time is proportional to L. Hence the center of mass diffusion constant D 
scales as 

D ~ LDo ~ L^'^-')'^ ^ L-^ (45) 

where the values z/ = 1/4 and 6 = 3/2 were used. The result D ~ L^^ is consistent with 



the simulations and intuitive arguments of Ref [^. It is remarkable that the diffusion 
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coefficient scaling is the same as in a linear reptating chain |T^, though the microscopic 



"walker" responsible for the motion has very different dynamics. In the case of a linear 
reptating chain, the exponents for the chain size and walker lifetime, v = 1/2 and 6 = 2, are 
different from the ring case, but they combine to give the same scaling for D. The longest 
relaxation time, in which the ring diffuses about one average ring size, scales as T^ax ~ L^^'^ 
(compare to the linear reptating chain result Tmax ^ Lr") ■ 



The relaxation behavior is very different from that of a linear reptating chain p5[. In 
the linear reptating chain, for most structures to relax, the head of the chain must reptate 
and free them, a process takes of the the order of the longest relaxation time Tmax ^ L^- In 
contrast, the hernia mechanism allows most structures to relax very quickly. 

5. Connection with Spontaneous Symmetry Breaking in ID: In this section, 
we discuss the models of gel-electrophoresis in the context of systems that display sponta- 
neous symmetry breaking (SSB) in one dimension (ID). It is well known that in thermal 
equilibrium, symmetry breaking and long range order can not appear in ID systems with 
short-range interactions. Recently, a simple non-equilibrium model that displays SSB in 
ID was presented |^^. The Rubinstein-Duke |^ reptation-tube model offers an additional 



interesting example of a ID system with short-range interactions which displays symmetry 
breaking. In this model, symmetry breaking has a particularly simple physical meaning. 
Consider a linear DNA chain with no tension or hernias. It is typically aligned with the field 
ip. There are two symmetric configurations: either segment number 1 is at the chain head, 
or segment number L is at the chain head. A finite chain eventually flips between these two 
configurations. The flipping time, however, grows very quickly with the chain size: in order 
to flip, the tail must reptate an order of L steps against the fleld, and the rest of chain must 
follow it, and thus the flipping time goes as ~ e°'^ with some constant a. The chain is 
therefore effectively stuck in one of the two configurations and symmetry is broken in the 
thermodynamic limit. 



In this example, as in the model of Ref P^, the mechanism permitting spontaneous 



symmetry breaking depends on the presence of boundaries (i.e. the existence of a chain 
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head and tail) and on a conservation law in the chain bulk (the number of + and — particles 
is conserved by the bulk dynamics when hernias are not permitted). The latter is seen to be 
important by allowing hernias, as in the present model (without tension which corresponds 
to long ranged interactions). Hernias amount to pair creation and annihilation in the bulk, 
and break the conservation law. Adding hernias also cancels the SSB: the chain is typically 
stuck in a symmetric hooked state (Sec 3a., Fig 5). 

We note that an interesting case of SSB in a ID system, with periodic boundary condi- 
tions and with no bulk conservation of the order parameter has been presented [^ . 
6. Discussion 

The role of tension and hernias in the dynamics of ring-shaped DNA gel-electrophoresis 
is studied. A microscopic model, which adds tension and hernias to the reptation model 
of Rubinstein and Duke, was proposed. The model is similar to react ion- diffusion models, 
with the addition of a new and interesting hierarchical pairing which represents the hernia 
fingers. 

We predict that the mobility of open-circular DNA should saturate much as in linear 
chains. Tension serves to stabilize a quasi-linear conformation of the ring. In this regime, 
there are many coiled DNA segments, screening the effects of tension. The screening of 
tension by coiled segments is essential for explaining the success of reptation models for 
linear DNA fragments [[I0| , p2| , which neglect both hernias and tension. It is suggested 



that tension and hernias have, to a certain extent, canceling effects. It is thus important to 
include both effects in models of hernia-mediated polymer dynamics. Some of the qualitative 
features of the motion of the ring are captured by a simple mean field theory. 

The present predictions are difficult to check in standard electrophoresis experiments in 
agarose gels in which the DNA rings become hooked on dangling ends in the gel ("hoop in 
stick" effect. Sec 3b). An experimental situation which eliminates the effects of dangling 
ends is offered by the microlitthographic arrays of posts introduced by Volkmuth and Austin 
P^ . These arrays were used to study electrophoresis of linear DNA, which could be observed 
by fluorescence microscopy. Since these 2D arrays have a floor and a ceiling attached to the 
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posts, a ring-shaped DNA can not become impaled by a post. Thus the hoop in stick effect 
is negated. The present work predicts that the mobihty of plasmids in such an obstacle 
array should saturate at a finite value, and not decrease exponentially with the ring size as 
in standard gels with dangling ends. 

In addition, agarose gels exhibit a very broad distribution of pore sizes. The present 
model pertains more closely to experiments on regular lattice, since the broad distribution 
of pore sizes in standard agarose gels might well affect the dynamical features and scaling 
laws discussed in this paper. 

The model is also applied to study the dynamics in zero field. The solution of the model 
in this case allows an exact foundation of some known scaling properties of ring polymers. 
The main effect of hernias on the static chain properties are to allow many undulations 
and wiggles per unit contour length. As a result, the ring adopts much more compact 
configurations (average size R ~ L^^^) than in the reptation case (where R ~ L^/^). The 
dynamics have been studied using a new approximate equation of motion for marked hernias. 
This equation seems to be a useful basis for analytical study of the effect of hernias on 
polymer dynamics. The lifetime of a hernia was found to scale as r ~ L^/^, which is shorter 
than expected from pure diffusion r ~ L^. This is due to an entropic bias of small hernias 
to shrink and be annihilated. The diffusion coefficient is found to scale as in the reptation 
case, D ~ L~^, though the behavior of the microscopic walker is different in the reptation 
and ring cases. The relaxation spectrum is very different, with most structures relaxing on 
short time-scales. These results were obtained for ring-shaped polymers, but should also 
be applicable to long linear chains in which hernia creation in the bulk is allowed. An 
exact solution of the steady state is also possible in this case, showing a crossover between 
reptation and hernia behavior as a function of polymer length. This solution, as well as a 
more detailed study of the zero-field dynamics, will be presented in a future publication. 

Finally, we note that the physically intuitive treatment of chain tension presented here 
may be useful also for describing different systems in which polymers in a melt or obstacle 
array move under external forces. 
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FIGURES 
FIG. 1. Configuration of a DNA chain (heavy line) in a gel, defined by a periodic lattice of gel 

pores (dotted lines). The DNA is divided into persistence length segments, numbered i = 1,2- ■ ■ L 

(in this case L = 8). The reptation tube (light line) is defined by all pores through which the 

DNA threads. The configuration is encoded using a + for segments that are stretched between two 

pores with the field direction, — for segments stretched against the field direction, and for coiled 

segments in the same pore. The displayed configuration is thus —0 + +00 — \- for i = 1 • • • 8. Move 

A corresponds to to a standard reptation move (inside the reptation tube) +0 -^ 0+. This move 

occurs with rate p. The reverse move OH > +0 is against the field direction and is given a smaller 

rate q. Move B corresponds to the formation of a hernia (leak through the reptation tube). It is 

represented by pair creation 00 -^ +u— • 

FIG. 2. A configuration of a ring-shaped DNA in the gel. The segments are numbered 
counter-clockwise. The configuration contains several branching fingers, and is encoded as shown 
by a string of +,— and O's along with their pairings into hernias . Ring-shaped DNA can move 
through the gel only by the annihilation and creation of hernias. 

FIG. 3. An example of a configuration of +,— and charges in the model, a) One of the 
allowed pairings into hernias, b) A forbidden pairing, where the hernias are not nested. Such a 
configuration can not be reached from an allowed configuration by the model dynamics. 

FIG. 4. Monte-Carlo simulation results of the model including hernias but neglecting tension, 
for a ring of size L = 52, with p = 0.56, q = 0.5 and c = a = 1. The mobility of the ring (center 
of mass velocity) is shown as a function of time (in sweeps, where one sweep equals L single-bond 
moves). The mobility is seen to have a spiked behavior, where the mobility is mostly zero with 
intermittent periods of motion. 
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FIG. 5. The ring configurations between one mobile burst and the next. The ring begins with 
a quasi-linear shape aligned with the field, which quickly develops an instability to secondary 
hernias and goes to a 2-armed hook. The hooks (in the absence of tension) are stuck for long 
times. Unhooking occurs by the retraction of one of the arms by fluctuations against the field, 
until a quasi-linear shape with a non-zero mobility is reached again. The spikes in the ring mobility 
in Fig. 3 are thus explained by the cycle between quasi-linear and hooked conformations. 

FIG. 6. Configuration of charges and hernia-pairings in the model that corresponds to a) 
Quasi-linear conformation, b) Hook with two equal arms. The quasi-linear conformation typi- 
cally also includes many sites, as well as small branching sub-hernias. 

FIG. 7. Mobility of a ring-shaped DNA fragment as a function of size. Monte-Carlo results for 
two field strengths, E = 2 (light lines) and E = 1 (bold lines) are shown with hernia creation and 
annihilation rates a = e~^ and c = 0.3. Units are such that Qb/kT = 1 so that the dimensionless 
field e = E. The mobility of the model including tension (full lines) decreases with the chain length 
for short chains, and then shows a saturation. For stronger fields, the asymptotic mobility is higher, 
and the chain length L* at which the mobility saturates is smaller. The results of the model without 
tension (dashed lines) are close to the results with tension for short chains (L < L*). However, 
without tension, an exponentially decreasing mobility is predicted for long chains, because of the 
formation of hooks. 

FIG. 8. Density of coiled segments, p, and mobility, /i, of ring-shaped DNA as a function of 
the ratio of hernia creation and annihilation rates c/a. The ring size is L = 40, the field strength 
IS E = 2 and the annihilation rate is held constant a = e~^. Shown are Monte-Carlo simulation 
results of the model including tension (o's), and the mean- field prediction (line). 
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